Splitting of the chiral critical point and realization of solitonic pion condensate driven 

by isospin density 
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We study the influence of the isospin asymmetry on the phase structure of strongly interacting 
quark matter near the tricritical point (TCP) using a generalized Ginzburg-Landau approach. The 
effect has proven to be so drastic, not only bringing about the shift of the location of TCP, but 
resulting in a rich fine structure at the vicinity of TCP. In particular, we find that an arbitrary small 
perturbation due to isospin density lifts the degeneracy of TCP making it split into four independent 
multicritical points. Accordingly, the homogeneous pion condensate and its solitonic counterpart 
come to occupy large domains in the Ginzburg-Landau coupling space. 

PACS numbers: 12.38.Mh, 21.65.Qr 



Introduction. — Quantum chromodynamics (QCD) at 
finite temperature and/or finite density is expected to 
exhibit a rich phase structure, and has been the sub- 
ject of extensive theoretical and experimental studies. In 
particular, several approaches to QCD with two light fla- 
vors suggest the existence of a critical end point (CEP) 
where the first order chiral phase transition turns into 
a crossover the CEP continues to a tricritical point 
(TCP) in the chiral limit where three lines of second- 
order phase transition meet. Despite many efforts based 
on the first principle calculations , the precise location 
of CEP in the phase diagram is still controversial. 

Recently based on a Ginzburg-Landau (GL) approach 
[H, and also in other effective models Q, it was shown 
that such TCP, if any, is actually replaced by a Lif- 
shitz point where three different forms of matter meet; 
i.e., a phase with broken chiral symmetry, a symmetric 
(Wigner) phase, and an inhomogeneous phase charac- 
terised by an additional translational symmetry break- 
ing. Such inhomogeneous chiral condensates can be 
viewed as a microscale ordered phase separation [f| @| . 

The effect of a finite quark mass on TCP is rather sim- 
ple; just to turn it into a CEP. Our focus here is the other 
important ingredient in realistic systems, i.e., the effect 
of an isospin asymmetry on TCP. Such a flavor symmetry 
breaking can be caused by a neutrality constraint which 
should be imposed in any bulk system to prevent the 
energy density from diverging. The isospin asymmetry 
effect on the thermodynamics is also important for the 
physics of heavy ion collisions with neutron rich nuclei at 
the Fermi energies, as it affects the isospin distillation in 
the nuclear liquid-gas phase transition [7|. 

The isospin asymmetry effects were well studied in the 
physics of color superconductivity at high density, and 
proved to lead to a rich variety of phases On the 
other hand, while there are a few model-based studies 
in the context of TCP [9|, there has been to present, to 
our knowledge, no systematic analysis based on the GL 
approach. Here we present for the first time GL analyses 
of the effects on phases at the vicinity of TCP. 



By comparison with other approaches, our GL frame- 
work has an advantage that it can give model indepen- 
dent predictions near TCP. Based on it, we take into 
consideration the charged pion condensate which may 
become favored by the inclusion of isospin chemical po- 
tential [10]. Moreover we incorporate the possibility of 
inhomogeneous condensate such as a chiral spiral |11| in 
which the chiral condensate is entangled with the pion 
condensate. Since we are interested in the response of 
TCP and its neighborhood against the isospin chemical 
potential (//i), our strategy is to take m as a perturbative 
parameter and expand the GL functional with respect to 
it. As a result we have many GL couplings, but it is pos- 
sible to derive universal relations among them; this will 
be performed based on the assumption that the quark 
loops are dominant near the TCP. This is justified if it is 
located at large fugacity region e M ^ T > 1. The analyses 
for the opposite case will be reported elsewhere. 

Our findings are; i) an arbitrary small perturbation 
due to isospin density not only brings about the shift 
of TCP but makes it split into four independent criti- 
cal points, and ii) both the homogeneous and inhomo- 
geneous pion condensates come to occupy large domains 
in the GL coupling space. Although it is beyond the 
scope of this Letter, these pion condensates near TCP 
may smoothly continue to pion condensates in nuclear 



matter at zero temperature [12| . or to those with some 
spatial structures [l3(; such meson condensates are of a 
renewed interest since their analogs may have a chance 
to be realized in ultracold atomic gases as was recently 
reported The fine phase structure near TCP may 

also have some phenomenological impacts on the physics 
of heavy ion collisions, and on the proto-neutron star 
cooling via the neutrino diffusion. 

Generalized Ginzburg-Landau approach at finite fj,i . — 
Let us start with writing the most general GL poten- 
tial for the chiral four vector </>(x) = (<x(x), 7r(x)) with 
ct(x) ~ —{tpi/}), 7r(x) ~ — (-0i7 5 x-0). We retain up to the 
sixth order in the order parameter and its spatial deriva- 
tive so as to allow for a minimal description of TCP. Then 
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the GL functional can be decomposed into three parts: 

Qgl[c(x), 7r(x)] = u>o + 5uj m + Suii, 

with wo being the chiral SU(2)r x SU(2)l invariant part, 
and Sujm (Slui) being the feedback from current quark 
mass (isospin density). The form of each part is rather 
stringently constrained by symmetry as we describe be- 
low. First, ujq can be set as 

-„ = f^ + f (0 2 ) 2 + ^(v0) 2 + ^W 

+ ^,V0) 2 + ^ 2 (V0) 2 -(0,V0) 2 ] 
+ ^(A0) 2 , 

where (</>, V(/>) denotes the inner product: cVer + ir ■ V-7T. 
Second, 8u>m — —ha breaks the chiral symmetry explic- 
itly: SU(2) L x SU(2) R SU(2) V . The GL coupling h is 
proportional to the quark mass m for light flavors. Lastly, 
the term Sojj represents the response to isospin density, 
which is our main focus here, 

r 02 _2 , 04 4 , 2 2\ 2 , Ate ^ 2 

6UI = —TT C + — TT C + —((/> - TT c )n c + — (V7T C ) , 

where ir c = (tti,^) is the charged pion doublet. This 
term breaks SU(2) L x SU(2) R down to U(1)i 3;V xU(1)i 3; a. 
Negative /?2, 04, 04b favor the homogeneous pion conden- 
sate |-7T C | 7^ while negative (3^ does inhomogeneous one. 
It is safe to neglect the quartic terms 04, /34b, 04c as long 
as we are only interested in the extreme vicinity of TCP 
where both |7r c |, \cf>\ are much smaller than the order of 
(ii. We here retain them because we are interested in the 
structure of phases where the order parameters become 
comparable with fii. To summarise the symmetry struc- 
ture, SU(2) L x SU(2) R is broken to SU(2) V due to 8lj m , 
and then further down to U(l)i 3: v via 8u>i. The residual 
U(l)i 3i v may be broken spontaneously via the formation 
of charged pion condensate. 

Ignoring the quark mass term Slom, we have eleven 
couplings {a 2 , a 4 , a ih , a 6 , a 6 ,6, a 6 , C! a 6jd , 02,134, 04b, 04c}- 
We can reduce the number of couplings significantly using 
the expansion about fii — as we demonstrate below. In 
order to make this, we assume that at the vicinity of TCP 
quark loops are dominant to gluonic ones, which may be 
justified provided that TCP is located at large fugacity 
region. The feedback of quark loops to the energy is 

Att = -^£-Tr(Sb£(x)f, 
V ^— ' n 

n=1 

where V denotes the volume of the unit cell of periodic 
condensate, and Tr should be taken over the Dirac, fla- 
vor and functional indices. So = diag (S u , Sd) is the 
bare quark propagator, and E(x) = ct(x)1 + i7 5 7t(x)t 1 
is the self-energy. Since the potential is symmetric under 
U(l)i 3; v x U(l)i 3; A that leaves a 2 + ir 2 and it 2 indepen- 
dently invariant, we set 7T2 = ttq = without any loss 




FIG. 1. The Feynman diagram for fiz, which is the leading 
driving force to the pion condensate when > 0. 



of generality to derive the GL couplings. The propa- 
gator in the momentum space is S(p) — diag^" 1 ,^" 1 ) 
with = (iuj m + fif,p) where the subscript / refers 
to the flavor index for u and d quarks, and U) m is the 
fermionic Matsubara frequency. We do not show all the 
results, but for example 02 and a 4 have the expressions; 

^=4TiV c E m , p i ^f £ ,a4 = 4TiV c E m , p ^. Let 
us start with the GL couplings a's which enter in ojq. 
From the expressions for a's, we can derive the universal 
relations ot4,b = a 4 , (ae.b, ae,c, ot6,d) = (5,3, l/2)ae 0). 
Although cjo is the chirally symmetric part, the couplings 
themselves can have dependence on /j, 2 , whose expansion 
should have the following general structure: 

fa 2 \ (I atf 0(rf)\ /a 2 0) \ 
a 4 = 1 bn 2 4°) , 

W \o o i J \ a f)j 

where a, b are just some constants and the subscript (0) 
denotes the quantity in the absence of isospin density. 
Explicit computation leads a = 0, b = 1. Similarly the 
expansion of Su>i is cast into the form: 

/ 2 \_ fen 2 0(f4) \ (ai 0) \ 
{04 [or 04 Hc) ]J \ drf [d b{c)l i\]) \af ) ' 

and we found c = —1/2, d = du c ) = — 1. When a® > 0, 
(3 2 = —fi 2 a i f > /2 acts as the main driving force to the 
formation of homogeneous pion condensate, whose dia- 
grammatic expression is shown in Fig. [1] Since we can 
use (> 0) to set the energy scale, we are left with only 
three parameters {a.^ , otf^ , A* 2 }■ We study the phase di- 
agram in the space of {a 2 \oif^} and how it is affected 
by /^j. It has in principle a unique map onto the QCD 
phase diagram at the vicinity of the critical point where 
the chiral condensate has a size less than or of order 

Just to avoid notational confusion, we suppress the 
subscript (0) for a's in the following. Moreover we scale 

— 1/2 

every quantity with energy dimension in the unit [a 6 ] . 
In particular we have ae = 1. Using the scaling argument 
as in H , it turns out that any critical condition is of the 
form /(a2/a 2 , /i 2 /|a 4 |) = where / is some function hav- 
ing two arguments. From this, we see that even when fij 
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is magnified as A/ii with A being an arbitrary real num- 
ber, the phase diagram stays the same once we redefine 
(a 2 j Q! 4) by (A 4 Q!2, A 2 a4). In particular, we conclude that 
the coordinate of any critical point in (0:2, a4)-plane, if 
any, should scale as a 2 oc /if, 0:4 oc /1 2 . 

Phase diagram with homogeneous states only. — Let us 
first discuss the consequence of nonzero /ii to homoge- 
neous condensates. We show the phase diagram in the 
lower panel of Fig. [2] Just for comparison, we also show 
in the upper panel how Su>m affects TCP and its neighbor- 
hood in the absence of fi\. In the latter case we actually 
find CEP at (2.28/i 4 / 5 , ■ 



15j . In contrast, we 



-2.25/i 2 / 5 ) 

find a drastic change of phase diagram due to nonzero 
fii in the former case. As expected, we find that pion 
condensate replaces a major part of the phase of chiral 
condensate in particular for 04 > 0. This is because 
P2 = —^1^4/2 < favors the pion condensation. In fact 
the second-order phase boundary between the Wigner 
and pion condensed phases can be derived from the con- 
dition of vanishing quadratic term: 



= 



dir 2 



a 2 - At 2 0:4/2 . 



Another notable point in the phase diagram is the shift 
of the location of TCP from (0, 0) to (0, -tf) labeled 
by TCP'. The shift itself is consistent with model-based 
studies [§] , and the location can be understood by noting 



n GL [<7,0] = ^ 2 
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where we see that a 2 = together with 014 + /.if = 
defines a new TCP. For the region of (a 2 < 0, en < 0), 
there is a competition between the chiral and pion con- 
densates, leading to a first order phase transition between 
two phases as depicted in the figure. The coordinate of 
characteristic points is listed in the Table HI 

Phase diagram with inhomogeneous states included. — 
We now address the question what is the impact of in- 
clusion of inhomogeneous phases. We restrict the crystal 
structure to one dimensional ones 0, We analyse 

three cases, (i) a(z) ^ 0, n(z) = 0, (ii) a(z) — 0, n(z) ^ 
and (iii) a(z) 7^ 0, ir(z) ^ 0. For the case (i) or (ii), we 
can solve analytically the Euler-Lagrange (EL) equation 
[H S E3 ; leading to a solitonic condensate 

o~{z) [or 7r(z)] = \fvh sn(kz; v) , 



TABLE I. Location of critical points; CEP, TCP' and P in 
Fig. H Q, R in Fig. El 
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FIG. 2. (color online). Phase structure for homogeneous 
condensates: Solid (dashed) line represents first (second) or- 
der transition, while double dashed line does the pseudo 
critical line for crossover. (Upper panel): The effect of h 
in the absence of u\. When h = , TCP is located at 
the origin; we have a second order critical line 02 = 
for q 4 > and a first order one q 2 = j^ot\ for Q4 < 0. 
Nonzero h turns TCP into a CEP whose location is shifted 
to (2.28/i 4/5 , -2.25h 2/5 ). (Lower panel): The effect of ui for 
h = 0. Asymptotic behaviours of the first order critical line 
separating the a phase from the n phase can be derived an- 

o 2 4 

alytically: a?2 — ► ~^ a 4 + 24 + C^Mi) as Q 4 — ► —00, while 

aL~. _l 1^,2 ao n . 0. 
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1 1 2 
2~Ct4 + as Q4 



with v denoting the elliptic modulus. In the case (iii), we 
can not solve the EL equation analytically. We take here 
the variational method instead. The variational state we 
consider here is one so called chiral spiral [ll| where a and 
tt are entangled: a(z) = mcos(qz), %{z) = msm(qz). 
This assumption is motivated by the study of the Gross- 
Neveu model [jjj . We compute the phase diagram by 
compering energies for the three cases. Resulting phase 
diagram is displayed in Fig. [3l We find no window for the 
chiral spiral, but notice a new fine structure appearing 
below TCP'. TCP' changes into a Lifshitz point from 
which a solitonic chiral condensate (<r(x)) expands be- 
tween the homogeneous a and Wigner phases; this is in- 
deed expected [3j. In this case, however, we notice that 
(r(x)-phase does not continue to the region of <C — Hi I 
a major part of solitonic a island is taken over by a 
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FIG. 3. (color online). Phase structure for fii 7^ 0. The 
second order critical line connecting TCP' to Q is 0.2 — § (0:4 + 
/if) 2 . The second order critical line linking TCP' to R is 02 = 
^(Q4+/i 2 ) 2 . The first order critical line between the <r(z) and 
7r(z) phases (a curve QR) requires a numerical computation. 



solitonic pion condensate 7r(x). This is caused by the 
term — (|ck4 | — ^i)(Vct) 2 in the potential where juf dis- 
favors <t(x). Since the quartic term affects thermody- 
namics via the square of its coefficient, ~ a| — 2//j|o!4|, 
the effect overwhelms the other ones at negative large 04. 
This makes 7r(x)-phase replace cr(x)-phase eventually. As 
a consequence there appear two additional multicritical 
points denoted by Q and R in the figure. The precise 
locations of these points are listed in the Table. Q] 

Let us briefly discuss the effect of quark mass. Since 
it breaks explicitly the chiral symmetry down to SU(2)y, 
the phase boundary between the Wigner and the a ^ 
phase would be smoothen out. We expect, however, 
that pionic phases will not be affected so dramatically 
since it is characterized by the spontaneous breaking of 
U(l)i 3; v and so is nothing to do with the explicit breaking 
of axial symmetry. The detailed analyses are now under 
investigation and will be reported elsewhere. 

In conclusion, we performed a systematic GL analysis 
on the effect of isospin asymmetry on TCP and phases of 
its neighborhood. By incorporating the effect of isospin 
density perturbatively at the leading order, we first de- 
rived the GL potential which works at the vicinity of TCP 
where the order parameters become comparable with /Ut. 
Based on it, we studied how the isospin asymmetry af- 
fects the phase structure. We found that it has several 
remarkable effects; it does not only cause a shift of the lo- 
cation of TCP, but also brings about the development of 
sizable region for homogeneous and inhomogeneous pion 
condensates. This in turn leads to the appearance of 
several new multicritical points. 

Let us finally make some speculative remarks about 
the phases found near TCP. First, although we found 
a first order phase transition between the solitonic c(x) 
and tt(x) phases, it may be replaced with a fine structure 



once we determine a suitable functional form for two con- 
densates solving the coupled EL equation; for instance, a 
phase with er(x) entangled with 7r(x) may show up. Sec- 
ond, the pion condensed phases near TCP may continue 
smoothly to those discussed in nuclear matter at low tem- 
perature [l2[. Lastly, multicritical points observed near 
TCP' may have some experimental signatures such as 
those discussed for CEP [19|. These clearly deserve fur- 
ther investigations. 
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